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The recently discovered GRB afterglow is believed to be described reasonably well by 
synchrotron emission from a slowing down relativistic shell that collides with an external 
medium. To compare theoretical models with afterglow observations we calculate here the 
broad band spectrum and corresponding light curve of synchrotron radiation from a power-law 
distribution of electrons in an expanding relativistic shock. Both the spectrum and the light 
curve consist of several power-law segments. The light curve is constructed under two limiting 
models for the hydrodynamical evolution of the shock: fully adiabatic and fully radiative. We 
compare the results with observations of 7-ray burst afterglows. 



Delayed emission in X-ray, optical and radio wavelengths has been recently seen in a few 7-ray bursts 
(Costa et al. 1997, Groot et al. 1997, Frail et. al. 1997). This so-called "afterglow" is described reasonably 
well as synchrotron emission from accelerated electrons when a spherical relativistic shell collides with an 
external medium (Waxman 1997a,b, Wijers, Rees & Meszaros 1997, Katz & Piran 1997). Previous analyses 
have described the spectrum and light curve only over a limited range of frequency and time. In this Letter 
we discuss the spectrum over a wide range of frequency and derive the shape of the light curve from very 
early to late times. We focus on the optical and X-ray emission, where synchrotron self absorption is not 
important, and we assume that the shell is ultra relativistic. We allow for both adiabatic and radiative 
hydrodynamical evolution. 

Consider a relativistic shock propagating through a uniform cold medium with particle density n. 
Behind the shock, the particle density and the energy density are given by Ajn and 47^nmpC^, respectively, 
where 7 is the Lorentz factor of the shocked fluid (Blandford & McKee 1976). We assume that electrons are 
accelerated in the shock to a power law distribution of Lorentz factor 7e, with a minimum Lorentz factor 
7m: N{'je)d'ye ot: 7~Pc?7e, 7e > 7m- To keep the energy of the electrons finite we take p > 2. We assume 
that a constant fraction Cg of the shock energy goes into the electrons. Then 



where the coefficient on the right corresponds to the standard choice, p = 2.5 (Sari, Narayan & Piran 1996). 
We also assume that the magnetic energy density behind the shock is a constant fraction of the shock 
energy. This gives a magnetic field strength 
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In this Letter we consider only synchrotron emission, and neglect inverse Compton scattering (which 
can be important when es > Ce, Sari et al. 1996). 



2. Synchrotron Spectrum of a Relativistic Shock 

A relativistic electron with Lorentz factor 7e ^ 1 in a magnetic field B emits synchrotron radiation. 
The radiation power and the characteristic frequency are given by (Rybicki & Lightman 1976) 

P(~le) = I'TClhl^, (3) 

where the factors of 7^ and 7 are introduced to transform the results from the frame of the shocked fluid to 
the frame of the observer. The spectral power, Pj, (power per unit frequency, ergHz"^ s~^), varies as u^/^ 
for ;y < z^(7e), and cuts off exponentially for v > z^(7e) (Rybicki & Lightman 1976). The peak power occurs 
at J^(7e), where it has the approximate value 

P « EhA - ih^^B (b) 

Note that Pv,max does not depend on 7e, whereas the position of the peak does. 

The above description of describes the emitted spectrum when the electron does not lose a significant 
fraction of its energy to radiation. This requires 7,. to bo loss than a critical value 7c, above which cooling 
by synchrotron radiation is significant. The critical electron Lorentz factor 7^ is given by the condition 
77cmec2 = P(7c)i, 

67rmeC Sm-e 1 

OT^B'^t IGeBcrTnipCtj^n 

where t refers to time in the frame of the observer. 



Consider now an electron with an initial Lorentz factor 7e > 7c- This electron cools down to 7c in the 
time t. As it cools, the frequency of the synchrotron emission varies as J/ oc 7^ while the electron energy 
varies as 7e. It then follows that the spectral power varies as over the frequency range i>c < < i^ile), 

where we have defined Pc = '^(ic)- The net spectrum of radiation from such an electron then consists of 
three segments, a low energy tail for v < Vc where goes as i/^/^, a power-law segment between Uc and 
i/(7e) where ~ and an exponential cutoff for u > z^(7e). The maximum emissivity occurs at Vc 

and is given by P„,max- 

As dcscribcid in the previous section, we are intercstcid in a power-law distribution of electrons. To 
calculate the net spectrum due to all the electrons we need to integrate over jg. There are now two different 
cases, depending on whether 7™ > 7c or 7™ < 7c. 

Let the total number of electrons be Ng. When 7^ > 7c, all the electrons cool down roughly to 7c and 
the flux at Vc is approximately N(.Pv,max- We call this the case of fast cooling. The flux at the observer, P^, 
is given by 

(j//j/c)l/3P. 

-Fl = •( {iy/iyc)~^^'^F^^max, Vm>y> Vc, (7) 

.max, V > Vrm 
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where Vm = z^(7m) and F^^^ax = -^efi/,max/47rZ)^ is the observed peak flux at distance D from the source. 

When 7c > 7^, only those electrons with 7e > 7^ can cool. We call this slow cooling^ because the 
electrons with 7^ ^ 7„j, which form the bulk of the population, do not cool within a time t. Integrating 
over the electron distribution gives 



Typical spectra corresponding to fast and slow cooling are shown in Figures la and lb. In addition to 
the various power-law regimes described above, self-absorption causes a steep cutoff of the spectrum, either 
as z/^ or z/^/^, at low frequencies (Katz 1994, Waxman 1997b, Katz and Piran 1997). For completeness, we 
show this regime in Figure 1, but we shall ignore it for the rest of this Letter since self-absorption does not 
affect either the optical or X-ray radiation in which we are interested. 



The instantaneous spectra described in the previous section do not depend on the hydrodynamical 
evolution of the shock. The only assumption made there is that the shock properties are fairly constant 
over a time scale comparable to the observation time t. The light curves at a given frequency, however, 
depend on the temporal evolution of various quantities, such as the break frequencies Vm and Vc and the 
peak power NePv,max- These depend, in turn, on how 7 and A^e scale as a function of t. 

We limit the discussion here to the case of a spherical shock of radius R{t) propagating into a 

constant surrounding density n. Clearly, the total number of swept-up c;lc;c;trons in the post-shock fluid 
is Nf, = AnR^n/i. We consider two extreme limits for the hydrodynamical evolution of the shock: either 
fully radiative or fully adiabatic. The radiative solution assumes that all the internal energy created in the 
shock is radiated. This requires two conditions to be satisfled: (1) the fraction of the energy going into the 
electrons must be large, i.e. eg 1, and (2) we must be in the regime of fast cooling, jc < 7m- If either of 
these conditions is not satisfied, i.e. if ^ 1 or 7,; ^ 7^, then we have adiabatic evolution. 

In the adiabatic case, the energy E of the spherical shock is constant and is given by 
E — IQnj^ nnipC^ /n (Blandford & McKee 1976, Sari 1997). In the radiative case, the energy 
varies as oc 7, where 7 ^ (R/L)-^. Here L = [17M/(167rmpn)]i/^ (Blandford & McKee 1976, Vietri 
1996, Katz & Piran 1997) is the radius at which the mass swept up from the external medium equals the 
initial mass M of the ejecta (We used 17/16 instead of 3/4 in order to be compatible with the adiabatic 
expression and to enable a smooth transition between the two); we write M in terms of the initial energy of 
the explosion via M = E/^q(?, where 70 is the initial Lorentz factor of the ejecta. 

In both the adiabatic and radiative cases, there is a simple relation connecting the radius of the shock 
R, the fluid Lorentz factor 7, and the observed time t: t — R/ct'y'^c, where the numerical value of ct lies 
between ~ 3 and ~ 7 depending on the details of the hydrodynamical evolution and the spectrum (Sari 
1997a, Waxman 1997c, Sari 1997b, Panaitescu & Meszaros 1997). For simplicity we use t = R/A'^'^c for all 
cases. We then have the following hydrodynamical evolution equations, 




(8) 



3. Hydrodynamical Evolution and Light Curves 




(9) 
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Fig. 1. — Synchrotron spectrum of a relativistic shock with a power-law distribution of electrons, (a) The 
case of fast cooling, which is expected at early times {t < to) in a 7-ray burst afterglow. The spectrum consists 
of four segments, identified as A, B, C, D. Self-absorption is important below Va- The frequencies, I'm, i^c, 
Va, decrease with time as indicated; the scalings above the arrows correspond to an adiabatic evolution, and 
the scalings below, in square brackets, to a fully radiative evolution, (b) The case of slow cooling, which is 
expected at late times {t > Iq). The evolution is always adiabatic. The four segments are identified as E, F, 
G, H. 
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, . ^ / (17^;/10247mmpc5^3)l/^ adiabatic, , . 

' ~ I (4c^/L)-3/^ radiative. ^ ^ 

Using these scalings and the results of the previous section, we can calculate the variation with time of all 
the relevant quantities. For the adiabatic case, we find 

l^c = 2.7 X 10l2£^3/2^-^l/2^_i^-l/2 

vra = 5.7 X lO^^.f .^£;,^/ Hz, 
F,,™a. = 1.1 X 10'eB^'^52«y'i?28' A^J, (H) 

where td is the time in days, £^52 = £'/10^^ ergs, ni is n in units of cm^'^ and D28 = I?/10^* cm. In the 
case of a fully radiative evolution, the results are 

- 1 V inl3^-3/2p-4/7 4/7 -13/14 -2/7 „ 
i^c - l.o X iU £5 72 ?^l irf Hz, 

- 1 9 V inl4^1/2 2p4/7 -4/7 -1/14 -12/7 TT 

P - ZL V in3^1/2p8/7 -8/7 5/14 „-2, -3/7 -[ „x 

where we have scaled the initial Lorentz factor of the ejecta by a factor of 100: 72 = 7o/100. 

The spectra presented in Figure 1 show the positions of Vc and for typical parameters. Note that, 
in both the adiabatic and radiative cases, Vc decreases more slowly with time than ^',„. Therefore, at 
sufficiently early times we have Vc < ^^m, i.e. fast cooling, while at late times we have Vc > i^m, i.e., slow 
cooling. The transition between the two occurs when Vc = Vm- This corresponds to the time 

^ r 210e^egi?52ni days, adiabatic, , , 

" \ 4:.Q€g^eJ^Et^2^^2^^^''^\^^ days, radiative. 

At t = to, the spectrum changes from fast cooling (Fig. la) to slow cooling (Fig. lb). In addition, if — ^ 1, 
the hydrodynamical evolution changes from radiative to adiabatic. However, if eg ^ 1, the evolution 
remains adiabatic throughout. 

Once we know how the break frequencies, Vc, Vnn and the peak flux F^^^nax vary with time, we can 
calculate the light curve. Consider a fixed frequency v = lO^^i^is Hz. From the first two equations in (11) 



and (|lj) we see that there are two critical times, tc and tm, when the break frequencies, Vc and Vm, cross 
the observed frequency v. 

_ J 7.3 X 10-^€g^E^2^ rij^ ^i^j^g^ days, adiabatic, 

= I 2.7 X 10-^6^2i/^i?5f 7|n^i3/V-^/^ days, radiative, ^ ' 

_ J 0.69e]/^ee^'^£^52'^i^]^^^'^ days, adiabatic, , . 

™ ~ 1 n on 7/24 7/6 ^^1/3 -1/3 -7/12 -1/24 , ,. l-'^^i 

[ 0.29e^ Ce £'52 72 ^15 "i days, radiative. 

It is easily seen that there are only two possible orderings of the three critical times, tm, to, namely 
^0 > tm > tc and to < tm < tc- Let us define the critical frequency, vq — Vcita) = i^mito), 

_ J 1.8 X 10^^e]f'^^e-^E^2^n^^^^ Hz, adiabatic, , . 

" \ 8.5 X l0^h-''^'"ec'/'E-2'/'^t^'n-''^'' Hz, radiative. ^ ^ 

When 1^ > vq, the ordering to > tm > tc applies and we refer to the corresponding light curve as the high 
frequency light curve. Similarly, when u < Uo, we have to < tm < tc, and we obtain the low frequency light 
curve. 
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Fig. 2. — Light curve due to synchrotron radiation from a spherical relativistic shock, ignoring the effect 
of self-absorption, (a) The high frequency case {ly > The light curve has four segments, separated by 
the critical times, tc, tm, to- The labels, B, C, D, H, indicate the correspondence with spectral segments in 
Fig. 1. The observed flux varies with time as indicated; the scalings within square brackets are for radiative 
evolution (which is restricted to t < to) and the other scalings are for adiabatic evolution, (b) The low 
frequency case {v < vo)- 
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Figure 2a shows a typical high frequency Hght curve. At early times the electrons cool fast and v is less 
than both Vra and v^- Ignoring self absorption, the situation corresponds to segment B in Figure 1, and the 
flux varies as Fy ^ F^,maxi'^/'^cy^^ ■ If the evolution is adiabatic, Fi^ ^aax is constant, and F^ ~ t^^^. In the 
radiative case, F^^max ^ and F^, ~ The scalings in the other segments, which correspond to C, 

D, H in Fig. 1, can be derived in a similar fashion and arc shown in Fig. 2a. 

Figure 2b shows the low frequency light curve, corresponding to v < vq. In this case, there are four 
phases in the light curve, corresponding to segments B, F, G and H. The time dependences of the flux are 
indicated on the plot for both the adiabatic and radiative cases. 

We conclude with two comments. First, during radiative evolution the energy in the shock decreases 
with time, the energy to be substituted in the radiative scalings is the initial energy. When a radiative 
shock switches to adiabatic evolution at time i = to, it is necessary to use the reduced energy to calculate 
the subsequent adiabatic evolution. The final energy -E/.52 which one should use in the adiabatic regime is 
related to the initial £'^,52 of the fireball by 

Ef,,, 0.022e^^/^e-3/5£;V5^-4/5^-2/5_ (^7) 

Second, if during the phase of fast cooling (t < tg) ee is somewhat less than unity, then only a fraction 
of the shock energy is lost to radiation. The scalings will be intermediate between the two limits of fully 
radiative and fully adiabatic discussed here. 

4. Discussion 

The main results of this Letter are summarized in Fig. 1 and Fig. 2, along with the scalings given in 
equations (|ll])-(|l|). 

It is well-known that the flux at the peak of the synchrotron spectrum is independent of time in the 
slow cooling limit for adiabatic hydrodynamic evolution (Meszaros & Rees 1996). We have shown in this 
Letter that the peak flux is constant even in the fast cooling limit if the evolution is adiabatic. The fast 
cooling stage has generally not been treated by other authors. We show that the position of the peak of the 
spectrum varies as Vc oc t~^^^ during fast cooling compared to Vm oc in slow cooling (Fig. 1), and this 

is one way of distinguishing between the two cases. We have also derived the scalings for a fully radiative 
evolution, where all the shock energy is radiated efficiently. This regime again has rarely been discussed in 
the literature. We find that the peak flux decreases with time as F^^max oc V'^l'' and the position of the 
peak varies as oc (These results differ from those given in Katz & Piran 1997, who considered the 

flux at Vra instead of the peak flux, which is at v^^ 

Even within the adiabatic case, we find that there are two possible slopes for the decaying part of the 
light curve. Writing the flux as Fy ~ i^", the two cases give a = 3p/4 — 3/4 and a = 3p/4 — 1/2. If the 
physics of particle acceleration in relativistic shocks is universal in the sense that the power law index p 
of the electron distribution is always the same, and if the evolution is adiabatic, then we expect always to 
observe one of these two values of a, which differ by 1/4. Indeed, some X-ray afterglows appear to decay 
with a = 1.4 while the optical and X-ray afterglows of GRB 970228 and GRB970508 had a = 1.2 (Yoshida 
et al. 1997, Sokolov et. al. 1997). The difference between the two values is consistent with 1/4. The 
corresponding value of p is ~ 2.5 — 2.6, which is a reasonable energy index for shock acceleration. If future 
observations of 7-ray burst afterglows always find decays with either a = 1.4 or a = 1.2, it will be a strong 
confirmation of the shock model and the adiabatic assumption. The characteristic values of a are different 
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for radiative evolution. 

In addition to the decay of the light curve with time, we can also consider the spectral index /3, defined 
by ^ . The two values of a given above for adiabatic evolution correspond io f3 — [p — l)/2 and 
(3 = p/2, respectively. Thus, the relation between a and (3 in an adiabatic fireball is either a = 3/3/2 
or a — 3/3/2 + 1/2. Previous studies (Rees & Meszaros 1996, Waxman 1997a) have considered the first 
possibility. However, we note that, for the standard choice of parameters, namely n ~ 1 cm"'^ (a standard 
interstellar medium), Ce, cb > 0.1 (rough equipartition of energy between electrons and magnetic fields), 
the second relation holds in both the optical and X-ray bands during much of the decay. Indeed, this 
relation is more compatible with detailed observations of GRB 970508 (Sokolov et. al. 1997). 

Finally, we note that in none of the cases considered does the flux rise more steeply than t^^^. This is 
a potential problem since GRB970508 displayed a sharp rise in the optical flux just before its peak at two 
days (Sokolov et. al. 1997). 
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